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Abstract 

The multiple scale expansion method is used to derive amplitude equations for 
a system with thermohaline convection in the neighborhood of Hopf and Taylor 
bifurcation points and at the double zero point of the dispersion relation. A complex 
Ginzburg-Landau equation, a Newell- Whitehead- type equation, and an equation of 
the if 4 type, respectively, were obtained. Analytic expressions for the coefficients 
of these equations and their various asymptotic forms are presented. In the case of 
Hopf bifurcation for low and high frequencies, the amplitude equation reduces to a 
perturbed nonlinear Schrodinger equation. In the high-frequency limit, structures 
of the type of "dark" solitons are characteristic of the examined physical system. 
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Introduction. 



In the 1980-1990s, a number of papers devoted to the formation of structures 
in the neighborhood of Hopf bifurcation points for systems with translational 
invariance along the horizontal appeared in the literature on double-diffusive 
convection. Oscillations in such systems can give rise to waves of various types 
(for example, standing, traveling, modulated, and random), which are conve- 
niently studied by constructing amplitude equations [lj. An amplitude equa- 
tion for a system with convection was first obtained by Newell and Whitehead 
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[2]. It describes two-dimensional thermal convection and has the form of a 
generalized Ginzburg-Landau equation. Coulett et al. [3j proposed a system 
of Ginzburg-Landau equations that describes traveling, double-diffusive waves 
propagating on both sides in a liquid strip which is infinite along the horizon- 
tal: 

(d t + sd x )A R = (c + ic x )A R + (c 2 + ic z )d 2 x A R - 

-(c 4 + ic 5 )\A R \ 2 A R - (c 6 + ic 7 )\A L \ 2 A R , 
(d t - sd x )A L = (c + ici)A L + (c 2 + ic 3 )d%A L - 

-(c 4 + ic b )\A L \ 2 A L - (c 6 + ic 7 )\A R \ 2 A L . 

The form of these equations is postulated from general considerations (such 
as considerations of symmetry); it is assumed that the coefficients in these 
equations should be derived by asymptotic methods from the partial equations 
describing a particular physical system. 

However, a thorough and well-founded derivation of amplitude equations for 
double-diffusive systems is not available in the literature there. In many pa- 
pers, the form of the coefficients in Eqs. ([I]) is not discussed. In some papers, 
these coefficients are obtained from various physical considerations. Thus, 
Cross [4], examining a system with convection for binary mixtures in the 
limit of low Hopf frequencies, set the coefficients ci, 03,05 and c 7 in Eqs. (JT|) 
equal to zero as a first, crude, approximation, motivating this by empirical 
data indicating an analogy between the case considered and the case of purely 
temperature convection. Clearly, such assumptions on the form of the coeffi- 
cients can be rigorously justified only in a rigorous mathematical derivation 
of amplitude equations. 

In papers on double diffusive convection of binary mixtures in bulk and porous 
media, the Hopf frequency turns out to be unity in the case of oscillatory con- 
vection. For thermohaline convection, it is reasonable to consider the asymp- 
totic behavior for the Hopf frequency tending to infinity. In this limit, the 
amplitude equation should become the nonlinear Schrodinger equations gov- 
erning internal waves in two-dimensional geometry. 

In the present paper, using the derivative expansion method, which is a version 
of the multiple-scale expansion method, we derive amplitude equations for 
double- diffusive waves in two-dimensional, horizontally infinite geometry in 
the neighborhood of the Hopf and Taylor bifurcation points and the double 
zero point of the dispersion relation. Idealized boundary conditions are used. 
In the case of Hopf bifurcation, the amplitude equation for waves propagating 
only in one direction is examined. Analytic expressions for the coefficients of 
these equations are obtained. Their various asymptotic forms are studied, and 
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asymptotic forms of the amplitude equations for various parameter values are 
discussed. 



1 Formulation of the Problem; Basic Equations. 

The initial equations describe two-dimensional thermohaline convection in a 
liquid layer of thickness h bounded by two infinite, plane, horizontal bound- 
aries. The liquid moves in a vertical plane and the motion is described by 
the stream function ijj(t,x,z). The horizontal x and vertical z space variables 
are used; the time is denoted by t. It is assumed that there are no distributed 
sources of heat and salt, and on the upper and lower boundaries of the regions, 
these quantities have constant values. Hence, the main distribution of temper- 
ature and salinity is linear along the vertical and does not depend on time. 
The variables 9(t,x,z) and £(t,x,z) describe variations in the temperature 
and salinity about this main distribution. There are two types of thermoha- 
line convection: the finger type, in which the warmer and more saline liquid 
is at the upper boundary of the regions, and the diffusive type, in which the 
temperature and salinity are greater at the lower boundary. In the present 
paper, we study the second type. 

The evolution equations in the Boussinesq approximation in dimensionless 
form are a system of nonlinear equations in first-order partial derivatives with 
respect to time that depend on four parameters: the Prandtl number a, the 
Lewis number r (0 < r < 1), and the temperature Rr and salinity Rg Rayleigh 
numbers [ofS] : 

(d t - a A) A^j + a (R s d x £ - R T d x 6) = -Jty, Aip), 

{d t -A)0-d x ij; = -J{ij;,e), (2) 

(«9 4 -rAK-a^ = -J(^0- 

Here we have introduced the Jacobian J(f,g) = d x fd z g — d x gd z f. The bound- 
ary conditions for the dependent variables are chosen to be zero, which implies 
that the temperature and salinity at the boundaries of the region are constant, 
the vortex vanishes at the boundaries, and the boundaries are impermeable: 

^ = 0^, = = f = O 2 = 0, 1. (3) 

In the literature, these boundary conditions are usually called free-slip condi- 
tions or simply free conditions since the horizontal velocity component at the 
boundary does not vanish. 

The space variables are made nondimensional with respect to the thickness of 
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the liquid layer h. As the time scale, we use the quantity to = h 2 /x, where x is 
the thermal diffusivity of the liquid. The vertical and horizontal components 
of the liquid- velocity field are defined by the formulas 



v z = -rd x il), v x = -—d z i>. 
n h 

The dimensional temperature T and salinity S are given by the relations 



T(t, x, z) = T_ + 5T [1 - z + 6(t, x, z)} , 
S{t, x, z) = S_ + 5S [1 - z + x, z)\ . 

Here 5T = T + — T_, 5S = S+—S-, and T + , T_ and S + , 5_ are the temperatures 
and salinities on the lower and upper boundaries of the region, respectively. 
The temperature and salinity Rayleigh numbers can be expressed in terms of 
the parameters of the problem: 



X^ 

Here g is the acceleration of gravity, v is the viscosity of the liquid, and a and 
7 are the temperature and salinity coefficient of cubic expansions. 



2 Dispersion Relation and Its Consequences. 



We consider a system of partial differential equations that is derived by lin- 
earization of the initial system (|2j) in the neighborhood of the trivial solution: 



(d t - a A) Ai> + cr(R s d x Z-R T d x t 
(dt-A)e-d x iJ> = 0, 
(d t -TA)£-a^ = 0. 



(4) 



These equations are solved subject to boundary conditions ([3]) by the method 
of separation of variables. We seek a solution in the form 



Ai exp(At — i/3x) + Ai exp(At + i/3x) 



sin (Ti z) 



(5) 



Here the bar denotes complex conjugation, ip is the vector of the basic de- 
pendent quantities cp = (i[),9,£), (3 is the horizontal wavenumber, and Ai = 
(qai, o>ti, a>si) is the amplitude vector. For qai we use the notation A = a^i- 
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Substitution of (j5J) into (J4]) gives a system of algebraic equations for the vari- 
ables a,Ai,aTi and asi- The condition for the existence of solutions of this 
system takes the form of an algebraic equation of the third order in A [5]: 



A 3 + (1 + r + a)k 2 \ 2 + [(r + a + ra) + a (r s - r T )} k A X + 

+cr (rs - rr T + t) k 6 = 0. 



(6) 



Here we introduced the wavenumber k 2 = ir 2 +f3 2 , and the normalized Rayleigh 
numbers r T = R T /R* and r s = Rs/R*, where R* = k 4 (k//3) 2 is the Rayleigh 
number, for which there is loss of stability of the steady state for purely 
temperature convection. 

Equation has three roots, two of which can be complex conjugate. In the 
physical system considered, loss of stability occurs when with change in the 
bifurcation parameters and rs, one or several roots pass through zero or 
gain a positive real part if they are complex. 

In the plane of the parameters Tt and rs (see Fig. [TJ), it is possible to dis- 



tinguish regions I and II, on whose boundary there is loss of stability. The 
boundary itself consists of two rectilinear segments. On segment 1, Taylor bi- 
furcation is observed when one of the roots of the dispersion relation passes 
through zero, which gives rise to steady drum-type convection. On segment 
2, Hopf bifurcation takes place when the real parts of two complex conjugate 
roots become positive. As a result, oscillatory convection occurs. The segments 
adjoin at the point C, at which the dispersion relation Q has a double root. 
At this point, the parameter values are defined by 
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Fig. 1. Plane of the parameters and rs- 



r T1 = 



a 1 — r 



1 r + a 



rsi = 




5 



The straight lines on which Taylor and Hopf bifurcations are observed, are 
given, respectively, by the equations 



r T = -r s + l, r T = 1 + -(1 + r + a) + — — r s . 
t a 1 + a 

The oscillation frequency of oscillatory convection is determined by the imag- 
inary part A and is expressed in terms of the reduced frequency Q as Im(A) = 
flk 2 , and Q is, in turn, calculated from the formula 

Q 2 = -t 2 + ^—^(rr s , A = ink 2 . 

1 + cr 

Below, the reduced frequency Q is called the Hopf frequency. 



3 Slow Variables and Expansion of the Solutions. 

We consider the equations of double- diffusive convection in the neighborhood 
of a certain bifurcation point, for which the temperature and salinity Rayleigh 
numbers are denoted by and R^, respectively In this case, the Rayleigh 
number is written as 

R T = R£ (l + e 2 V ) , R s = R* s (l + e 2 Vs ) . 

The values of r\ and r\s are of the order of unity, and the small parameter e 
shows how far from the bifurcation point the examined system is located. To 
derive amplitude equations, we use the derivative expansion method of [7118] . 
We introduce the slow variables 

T x = et, T 2 = e 2 t, X x = ex. 

Next, into the basic equations ([25, we introduce the extended derivative by 
the rules 

d t ^d t + ed Tl + e 2 d T2 , d x -> d x + ed Xl ■ (7) 
The dependent variables are represented as series in the small parameter: 

cp = £ e n <p n (x, t, X u T u T 2 ) + 0(e 4 ). 

n=l 
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Substituting these expressions into Eqs. (j2J) with derivatives extended accord- 
ing to ([7]) and grouping terms with the same power of e, we obtain 



0(e) : LVi = 0, 

0(e 2 ) : L*cp 2 = -(L x d Tl ~ L 2 d Xl )<Pi - Mi^, 
0(e 3 ) : L*<p 3 = -(L 1( 9 Tl - L 2 <9 X > 2 - 

-(Lid T2 + L 3 d Xl + L A d Xl d Tl + L 5 )cp 1 - M 2 ((p 1 , tp 2 ). 

Here the operators L\,L 3 and L4 have diagonal form: diagLi = (A, 1,1). 
diagL 3 = (d t — 2aA — Acrd^, —1, —1), diagL 4 = (2<9 X ,0,0); the operators L* 
and L 2 can be written as 



L* 



V 



(dt - cjA)A -aR* T d x aR* s d x 
-d x (d t - A) 
-d x (d t - r A) 



2aA)d x aR* T -aR* s ^ 



\ 



-2(d t 



2d x 
2rd x 



In the operator L 5 , only the upper row is different from zero: (0, — aR^T]d x , 
(jRgT]sd x ). The vectors Mj = (M^, M?i, M$i) with nonlinear terms have the 
components 



M Al = J(fa, AVi), M T1 = J(ViA), Msi = J(Vi,a), 
M A2 = J(</> 2 , A^i) + J(Vi, A^ 2 ) + J(^i, 2d x d x M + 

+d z Aip 1 d Xl ipi - d z ipid Xl Ail) u 
M T2 = J(V>i, 2 ) + 0i) + dAdx^i ~ d,ij> x dxA, 
M S2 = J(V>i,6) + J(V>2,£i) + ^i^i - dtfadx^. 

The systems obtained can be written in general form: 



Here the functions Q- consist of terms that are in resonance with the left side 
of the equations. The condition of the absence of secular terms in solutions 
of similar systems of equations is known (see [7f8] ) to reduce to the require- 
ment that the functions Qj and the solutions of the conjugate homogeneous 
equation (L*)*ip* = be orthogonal. We now derive the relation to which 
the condition of the absence of secular terms reduces in this case and which 
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will be used below to derive amplitude equations. Let us consider the inhomo- 
geneous system of algebraic equations that is obtained from (J4]) by choosing 
the single-mode form fl^D and substituting functions Qj = (Qau Qn, Qsi) 
such that Qj = Qj exp(At — i(3x) + Qj exp(At + i(3x) into the right side of the 
homogeneous system: 



(A + ak 2 )(-k 2 )a Ai + aR* T i(3a T i - aR* s i(3a Si = Q Ai , 

(X + k 2 )a Ti + i(3a Ai = Q Ti , (8) 

(A + rk 2 )a Si + if3a Ai = Q Si . 

The solvability condition for this system of equations is formulated as the 
requirement that the right side be orthogonal to the solution of the conjugate 
homogeneous system [9] (1, — z/5aR^/(A + k 2 ), i(3crR* s /(X + rk 2 )) and reduces 
to the equation 

(A + k 2 )k 6 ar* s Q Si - (A + rk 2 )k 6 ar* T Q Ti - 

-(X + k 2 )(X + rk 2 )if3Q Al = 0. (9) 

For A = 0, this relation is simplified: 



1 i(3 

-r* s Qsi - r* T Q T i tjQai 

t ak 4 



4 Derivation of Amplitude Equations. 



Let us assume that the solution of the equations for ip x has the form ([5]) 
and the amplitude of this solution now depends on the slow variables: A = 
A(Ti, Xi, T 2 ). Substituting this solution into the equations for cp 2 we obtain a 
system of equations of the form (IE]) in which the functions Q2 are written as 



For this algebraic system to be solvable, it is necessary that condition Q) be 
satisfied. At different bifurcation points, this condition is formulated as dif- 
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ferent equations. We consider successively the equations obtained from the 
solvability condition of the indicated system at the bifurcation points charac- 
teristic of the physical system considered. 

In the last relations, we substitute the value of A at the Hopf bifurcation point 
A = iVlk 2 and set k 2 / f3 2 = 3 and (3 = tt/\^2, which is valid for the oscillation 
mode that is the first to lose stability [15] . In addition, we take into account 
the relations 

1 a + r /„„ \ 1 a+ 1 



fi 2 + l), r* = --^±{tt 2 + T 2 



1 a 1 -t V J ' a al 

Then, Eq. ([9]) can be written as 



d Tl A + V27rnd Xl A = 

and solved in general form by introducing the new slow variable X = X\ — 
\/27rOTi. If we assume that the amplitude A(X,T 2 ) depends on Xi and Ti 
only via X, this equation becomes an identity. 

In the other cases, where we consider the system at the Taylor bifurcation 
point or at the double zero point, the solvability condition (jUJ) has the form 



-(1 - r) (4 - -p—j J d Tl A + 2 i p^-~3j d Xl A = 0. (10) 

If in this equation, as above, k 2 /f3 2 = 3, i.e., the least stable oscillation mode 
is considered, then 3t x A = holds for the case of Taylor bifurcation. In the 
case of the double zero point, Eq. ( fTUj) is satisfied identically. 



5 Amplitude Equation at the Hopf Bifurcation Point. 



We now write the solution for ip 2 with the wavenumber for which there is loss 
of stability of the steady state: 



<£2 



A 2 ex\){iVLk 2 t — if3x) + A 2 exp(— iVtk 2 t + i(5x) 
+B 2 sin(27T2). 



sin(7T2;) + 



Here A 2 = (a^, ay 2 , a S2 ) and B 2 = (0,6^2,^2) are vectors that depend on 
the slow variables. The components of these vectors have the values 
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J T2 



J_JAf_ 

6?r i + n 2 ' 



>S2 



1 \A\ 2 



2 1 /_ . 3i7T \ 

as2 = 9^7T^ \ dxA -T2 aA2 )- 



Using the solutions given above, we formulate a system of equations from 
which it is possible to find ip 3 . This system of equations, as the system for ip 2 , 
has the form (jSJ). Then, we can write the functions Q3 as follows, retaining in 
them only terms with A{X,Tz): 



Q M = ^tt 2 {d T2 A - i(4ifi + 7a)d 2 x A- 



3tt 2 



Qt3 
QS3 



2(1 -r) 
iV2 1 



[(o- + l)(r - zfi)»75 - (a + r)(l + ^77] A , 



3tt 

iV2 1 



1 - ■ ^ 1 A\A? 



37r r + 



7T 



a T2 ^--(2zO + 5r)^ + ^- 



A\Af 



Condition (jUj) for system (jSj) has the form 



(a + l)(r - ^)Q 53 - (o- + r)(l - zfi)Q T3 - (1 - r)^Q A3 = 0. 

After transformations, we find that the amplitude A(X, T2) should satisfy the 
complex Ginzburg-Landau equation 



d T2 A = a x A + (3 X A\A\ 2 + 7l d|A 



The coefficients in this equation are defined by the formulas 



Oil 



Pi 



3i7r 2 [^(a + l)(n 2 + T 2 )(itt + 1) - V (a + r)(fi 2 + 1)(^ + r)} 



111 

8Q 



71 = iQ + 2 



4fi[ifi + (H-«7 + r)](l-r) 
a + ffr + r)fi — icrr 



n[ifi + (1 + a + r)] ' 
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6 Equation in the Form of a Perturbed Nonlinear Schrodinger 
Equation. 



For further investigation, the equation obtained can be brought to a more 
convenient form. We set rjs = 0. This implies that the behavior of the system 
can be controlled by changing the temperature gradient in the layer, while 
the salinity gradient remains constant and equal to the critical value. The 
coefficient an (iai/rj = acR + iaci) is eliminated from the equation by changing 
the dependent variable by the formula A = A'exp(— iaj^T-i). Equation (fTTj) 
then becomes 



id T2 A' + ~i R d\A' - (3 R A'\A'\ 2 = ia^A' + ijjd^A'. (12) 

Here 



ocr 



li 



3 9 cr + r 

-7T 

4 1 



Q 2 + l 



-7T 



t VI 2 + (1 + r + a) 2 v 
<7 + t)(<7 + 1) n 2 + i 



t(1 



n 



7T 

8fi' 



1 - 



T 



tt-2 



n 2 + (l + r + a) 2 ' 

a + ar + r)fi 2 + ar{\ + r + a) 



n\n 2 



T + cr) 



(o- + r)(l + r + a + ra) 
n 2 + (l + r + a) 2 ' 



Thus, the amplitude equation becomes a nonlinear Schrodinger equation (NSE) 
with perturbation on the right side. We note that the coefficient values in this 
equation do not coincide with those given by Bretherton and Spiegel [10] , who 
studied an equation of the type (TTTT) obtained by the method of expansion of 
a linear dispersion relation in the neighborhood of a critical wavenumber. 

It is known [llj that for a dissipatively perturbed NSE, perturbation changes 
the form of solutions only slightly. As a result, as the perturbing terms tend to 
zero, these solutions become solutions of an NSE with no perturbed right side, 
which can be solved in general form by the method of the inverse scattering 
problem [7j. If an NSE has soliton solutions (both envelope solitons and soli- 
tons above a field of finite density), Eq. (|12p has solutions of the same form 
with rather small perturbing terms. The type of NSE is determined by the 
sign of the second derivative. In this case, the sign of the coefficient 7^ varies 
with change in Q from zero to infinity. Hence, in the problem considered, two 
types of NSE are possible: 
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-for n 







IR 



2ot 
1 + r + O 



-i 



2to 



t + o) 



li = 2 + 



1 + r + a 



TO — 1 



1 + t + a) 3 

TO 



1 + r 



o 



+ o(n 2 ) ; 



- for — > oo 



7iJ = - 2(r + o + ra)^ 1 + 0(fT 3 ), 

7/ = 2(r + a)(l + r + a + ra)^ 2 + 0(fT 4 ). 

In the limit Q — 0, the coefficient 7 r becomes infinity and Eq. f[T2"j) loses mean- 
ing. This limiting case corresponds to the double zero point of the dispersion 
relation. The amplitude equation in the ^-neighborhood of this point will be 
deduced below. As Q increases from zero to infinity, 7r changes sign, whereas 
7/ decreases monotonically, remaining always positive. The frequency Q at 
which 7/j vanishes is determined from the formula 



^ = ^(l + a 2 + r 2 ) 

For rather large o or small r, this formula has the asymptotic form SIq « 
2to(1 + t + o")/(l + r 2 + cr 2 ). In the case, when o = 7, r = 1/81. we have 
f^o ~ 0.1663778362. For the other values of o and r see Fig. [61 In the limit of 
the infinite o it is true £Iq = y/2r. 



1 + 



Iot{1 + t + o) 



l + 2 + T 



2\2 



7 Transformation to a Nonlinear Schrodinger Equation. 

We consider two cases where the amplitude equation derived above becomes 
an NSE. Using the substitution 

A = yjaj/ '/3 R exjp[-i(a R + a I p 2 )T 2 ]F(a I T 2 , \Jai/i R X), 
where p is a positive constant, we bring Eq. f|T2|) to the form 

iF T + F xx - F(\F\ 2 - p 2 ) = i V F + ipF xx , (13) 



14 



0.2 








2 



4 



Fig. 6. Frequency 0,q(t, a) at which vanishes, for r = 1/81 (solid line), r = 1/10 
(dashed line), r = 1/4 (dotted line) and <r in the range — 5. 

where /i = ji/jr (Fig. [7]). Here and below, the subscripts T and X denote 
partial derivatives with respect to the slow time T2 and the X coordinate, 
respectively. The coefficient p tends to zero with increase in Q according to 
the asymptotic relation /i^2(t + o~)(1+t + o~ + ra)Q~ 3 . In addition, in the 
immediate vicinity of the Hopf bifurcation point (in the e 3 neighborhood), the 
first term on the right side of Eq. (fl~3l) can be eliminated. The second term 
can also be ignored if the frequency fl is sufficiently high. As a result, Eq. ([121 
becomes the NSE 



iF T + F xx - F{\F\ 2 - p 2 ) =0. 

This equation has solutions that are known as solitons of finite density or 
"dark" solitons [12]: 



The parameters ( and Xq characterize the width and initial position of the 
soliton, respectively. 

Thus, the present investigation shows that for the physical system considered, 



F = p 



exp«) + exp$ 2 2 f sin 2 (C/2) 

, Lr = 1 k 

1 + exp $ \ cosh $ 

$ = -pTsinC±(X-X )y27sin(C/2). 



) 



(14) 
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100 i 




-100 J 

Fig. 7. Coefficient /i(O) in Eq. ([13]), a = 7, r = 1/81. 

0.05 0.1 0.15 




Fig. 8. Coefficient /j,(Q) in Eq. CC| for a small ft, a = 7, r = 1/81. 

along with other solutions, there can be solutions of the type of "dark" soli- 
tons, and this is true in the limit of high Hopf frequencies. Apparently, double- 
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diffusive convection at high Hopf frequencies can occur in ocean systems. An 
example of these systems is a so-called thermohaline staircase [13]. Inversions 
of a thermohaline staircase often have stratification parameters, which corre- 
spond to the beginning of diffusive convection, and the Hopf frequency Q is 
of the order of 10 3 -10 5 . 

When the Hopf frequency tends to zero, Eq. (|12p takes a different asymptotic 
form. In this case, we set 



Thus, /i — > as Q — > (Fig. [8]). As in the previous case, the first term on 
the right side of the equation can be eliminated by assuming that the system 
is in the immediate vicinity (in the e 3 neighborhood) of the Hopf bifurcation 
point. Then, again, Eq. (1121) takes the form of an NSE: 

iF T = F XX + F\F\ 2 . 

This equation has well-known solutions in the form of envelope solitons. 

It is interesting that localized wave packets, with which soliton solution can 
be compared, were observed in experiments on convection of binary mixtures 
at rather low Hopf frequencies (see, for example, [TlfTo] ). 



8 Equations at the Taylor Bifurcation Points and Double Zero 



We consider the case of Taylor bifurcation or bifurcation to steady roll-type 
convection. On the straight line on which this bifurcation occurs, the dispersion 
relation has a first-order root. For terms of the order of 0(e 2 ), the equation has 
the form d^A = 0, i.e., the amplitude does not depend on the slow variable T\. 
For terms of the order of 0(e 3 ) of the functions Q3, we obtain the expressions 



A = )fa I Jfaexp(-ia R T 2 )F(a I T2, -^fo^pj^X). 



Then, 



iF T - F xx - F\F\ 2 = irjF + %\iF X x 



where /1 has the following low-frequency asymptotic form: 




Point. 
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QA3 — 


Q 

y 4r 
-o-TT [r, 

4 


Qt3 = 


in / 
6^ V 


Q S3 = 


in 


6t 2 x/2 



fas - 77) - ?7s]^ + ^|«9t 2 A - ^c^A 

Substituting these formulas into the compatibility condition, we have the am- 
plitude equation 



d T2 A = a 3 A - (3 3 A\A\ 2 - l3 d 2 Xi A, (15) 
where 



3 2 r^jris -y)-ys R n 2 r* T (l - r 2 

"3 = ^ T —TTi \ 7i — I 7 — \i P3 " 



r* T {\ - t) - (1 + r/cr) ' ^ 4r rf(l - r) - (1 + r/tr) ' 

73 



4r 



^ (1 - t) - (1 + r/a) ■ 

This equation is similar in form to the equation derived in [2] and reduces to 
it if a salinity gradient is absent. 

We consider the e 2 neighborhood of the point of intersection of the straight 
lines on which Hopf and Taylor bifurcations are observed. At this point, the 
dispersion relation has a second-order root (Takens-Bogdanov bifurcation). As 
noted above, for the case of the most unstable convective mode, the equation 
obtained for terms of the order of 0(e 2 ) is satisfied identically. Therefore, it is 
not necessary to use the variable TI or to introduce other slow variables. For 
terms of the order of 0(e 3 ) of the functions Q 3 , we obtain the expressions 



Qm = 4(137) K CT + Vis - t 1 + °I T )V\A - ~jf Xl d Tl A - —d 2 Xl A, 

After substitution of these expressions into the condition of the absence of 
secular terms, we obtain the equation 



d 2 1 A-c 2 d 2 Xi A = a 2 A + (3 2 A\A\ 2 , (16) 
where 
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6it 2 <jt 





(l + a/T)ri-(l + a)ris 
(l-r)(l+r + <r) 



l + r + a 



Equations of this type are known as y? 4 -equations, and they cannot be inte- 
grated accurately by the method of the inverse scattering problem [7j. Some 
papers [T6fl7|ll8j consider amplitude equations at the double point for the con- 
vection of binary mixtures. According to |19j . the results obtained for thermo- 
haline convection are extended to the case of convection of binary mixtures, 
where it is necessary to allow for the thermodiffusion effect. Therefore, for the 
last case, all the equations at bifurcation points derived in the present paper 
are valid with the parameters of the problem converted accordingly (Prandtl, 
Lewis, and Rayleigh numbers). Knobloch [18] obtained an amplitude equation 
at the double zero point that has the form d\ A = CiA + C2A\A\ 2 in the main 
order in e {C\ and C 2 are constants). Equation (TIB]) can be regarded as its 
extension to the case of spatial modulations. Brand et al. [T6] gives another 
amplitude equation at the double zero point, which includes a term with a 
third derivative of the form dtd 2 A. Therefore, it differs from the equations 
derived by the multiple- scale expansion method used in the present paper. 



9 Conclusions. 

1. The derivative expansion method is used to derive amplitude equations for a 
system with thermohaline convection in the neighborhood of the main bifurca- 
tion points characteristic of this system. In particular, within the framework 
of a unified approach, we obtained the complex Ginzburg-Landau equation 
(fTTT) in the case of Hopf bifurcation, the Newell- Whitehead equation ffl5|) in 
the case of Taylor bifurcation, and Eq. ffl6|) of the <^ 4 type in the neighborhood 
of the double zero point of the dispersion relation. 

2. Analytic expressions for the coefficients of the equations considered are 
given. For the equation in the neighborhood of the Hopf bifurcation points, 
the formulas specifying its coefficients refine the previous results of |10j . For 
the other two equations, such formulas, as far as we know, have not been 
previously reported in the literature. 

3. It is shown that, for low and high frequencies, the amplitude equation in 
the neighborhood of the Hopf bifurcation points reduces to the perturbed 
nonlinear Schrodinger equations f[T2l with characteristic solutions in the form 
of envelope solitons. In the high-frequency limit, the type of "dark" solitons 
(Ti~4"l) are characteristic of the examined physical system. 

4. The equation of the type of (p 4 derived at the double zero point of the 
dispersion relation can be regarded as an extension of the equation obtained 
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in [18] to the case of slow spatial modulations of the amplitude. 
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